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Abstract 

The recent discovery that some of the coefficients of the viscosity tensor 
are negative is shown to invahdate the hydrodynamic approach to the vortex 
hquid phase of a type-II superconductor. A satisfactory theory requires 
retention of all the spatial gradients of the velocities or electric fields and 
not just the first derivatives, as assumed in a hydrodynamic theory. We 
illustrate such a procedure by using time-dependent Ginzburg-Landau theory 
to determine the electric field distribution near a single "twin-plane boundary" 
due to a current passing through the boundary. 
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There exists an extensively developed phenomenological hydrodynamic theory |T|-§[ for 
superconductors in the vortex liquid regime. In the recent work of Mou et al. ||^, the 
coefficients of the viscosity tensor were calculated from a more fundamental theory, in this 
case Ginzburg-Landau theory, and some were found to be negative. We show explicitly that 
this invalidates the hydrodynamic approach. We find that a satisfactory theory requires 
retention of all the spatial gradients of the velocities and not just the first derivatives, 
as assumed in a hydrodynamic theory. We illustrate our procedure by calculating the 
electric field distribution near a single "twin-plane boundary" when a constant current is 
being passed through that boundary. The electric field decays to its bulk value; whereas, 
the hydrodynamic theory with negative viscosity coefficients (as calculated) has unphysical 
oscillatory behavior of the electric field as a function of distance from the boundary. 

Because the motion of extended fiux lines plays a vital role, the transport properties 
of type-II superconductors near and below the Hc2{T) line are expected to be nonlocal. 
For instance, the current j produced by a small electric field E, is given in terms of the 
conductivity tensor, as: 

J,{r) = J a^,(r,r') E,{r') dr' (1) 

where 0"^,^ is nonlocal, i.e. cr^;/(r, r') 7^ cr^j^(5(r — r'). Assuming translational invariance, one 
has then in momentum space: 

j^(k) = a,M EM, (2) 

where nonlocality now implies o"^;y(k) 7^ Const. The nonlocality of the conductivity has 
been recently investigated experimentally. [§] 

The "hydrodynamic" approach assumes that an adequate description of the long- 
time, large-distance behavior can be derived from a truncation of the small-wave-vector, 
small-frequency expansion of the conductivity to the lowest non-trivial terms. For instance, 
the dc hydrodynamic conductivity tensor would be taken in the hydrodynamic approach to 
be: 
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cT^^(k,u; = 0) = cr^!.(0) + S^a^y kakp, (3) 

(in the notation of Mou et al. |]^). Moreover, hydrodynamics includes a tacit assumption 
that this truncation resuhs in a stable theory; Szzzz^ for example, must be positive. (We 
take the external field B to be in the z direction.) 

However, recent calculations of Mou et al. of the S"s from Ginzburg-Landau theory 
found some of them to have signs opposite of those required for stability. We extend their 
calculations and show that stability is regained if the full expansion in powers of the wave 
vector is considered. Then we proceed to calculate the zz component of the conductivity 
in the presence of a boundary and with it, we calculate the electric field for the case of a 
constant current jz being passed through a "twin-plane boundary" (at z = 0). The results 
are compared with those from the hydrodynamic theory. This example also illustrates 
the necessity for stability of including the normal component of the conductivity (assumed 
to be local). 

Our starting point is the time-dependent Ginzburg-Landau equation: 



with noise correlations: 



^(r,t) + ^ |^(r,t)|2^(r,t) = r/(r,t), (4) 



2k T 

(r/*(r, t) r/(r', t')) = ^ 5(r - r') 6{t - t') (5) 

and with mi 2 = and = rric. We use the symmetric gauge 

. —By , . Bx 

= and Ay = — . (6) 

As the parameter F sets the time scale, time will here after be measured in units of F; h 
will be set to one. The following calculations are based on the linearized version of eq. (^, 
and so strictly speaking above the Hc2{T) line; however, making the Hartree approximation: 

a^ff = a + b (IV^P) (7) 

can take one below the Hc2{T) line. 



First, let us sketch out the derivation of the bulk conductivity, returning to the case with 
a boundary afterward. The linear response to the field can be calculated via the Kubo 
formula which relates the conductivity to the current fluctuations as follows: 
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d(r - r') / dit - t') e^k.(r-r')-^..(^-0 ^ t) J.(r', t'))c • (8) 



2kBT 

Here, we will be interested in the non-local nature of the dc conductivity: a^^^(k,uj = 0). 
(The frequency dependence of the uniform conductivity cr^^J(k = 0,0;) was studied in ref. 
0.) 

In order to calculate the current fluctuations, we will need the Green's function and 
correlation function. Instead of the usual eigenfunction expansion, we take the result for 
the Green's function from the path- integral approach [pJj] : 



G(r,t' + r;r',t') = ( 



me 1 r rric ,2 

1 -^TTT- — 1-T exp <^ — ar (z — z I 

V327r3£V/ sinh(c<;or/2) 2r ^ ^ 



exp 



coth(co'or/2) 
4£2 



{x-x'Y + {y-y'Y\ + —{xy'-yx') , (9) 



where cOo = e*B/mab is the cyclotron frequency, i = {e*B) is the magnetic length, and 
r > 0. 

From the definition of a Green's function: 



^(r,t) = Jdt' J dr' G{r,t]r',t') r]{r',t') 



(10) 



which one can then use to calculate the correlation function {ilj*{r',t') ip{r,t)). Using eq. 
(^) for the noise correlations, it follows after some algebra that: 



roo 

{ij*{r',t') iP{r,t)) = 2kBT ds G{r,r';\t - t'\ + 2s). 

Jo 
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The Kubo formula for the conductivity involves a product of correlation functions with 
some (covariant) derivatives acting on it, those derivatives coming from the usual definition 
of the current density: 
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_d_ 



e*AAiP + \i 



_d_ 

dr., 



(12) 



The next steps of the calculation are relatively straightforward: (i) insert eqs. (|Tl|) and 
(p!^) into the Kubo formula (ii) carry out the required derivatives; and (iii) perform the 
spatial integrations in eq. which are all Gaussian integrals. Three integrations over 
time variables (one from the Kubo formula and one from each correlation function) remain. 
After some rearrangement of these time integrations, we find: 

f- (^'O) = F^T^ ^/ T7 \ ^^^/ du C^,{T,u,v,k) 

STT-^'^mnm^i^ Jo smh(ct;or) Jo Jo 



1 _ £^[cosh(.;or)- coshH^or)] , , _ r{l - uW) | 
1 2 smh(u;o7") ^ 4mc J 



where 



^ 2 cosh(ut>co'oT) sinh {uvujqt) ^ [cosh(co'oT) — cosh(Mt>co'oT)] ^2 
£2 sinh(co'oT) sinh^(a;or) sinh^(a;or) ^ 

r 

sinh^(Mt>ct;or) — [cosh(ci;oT) — cosh(Mwc<JoT)]^ 
° sinli2(^„T) ' " 

= fefc. (14) 

smh(C(JoTj 

The {a; ^ y\ symmetry can be used to obtain the other C^v's,. 

Expanding eq. ( jlSD to quadratic order in k and performing the u and v integrations 
reproduces the results of Mou et al. including results such as S^xxx < and Szzzz < 0. 
However, one can see from eqs. (^) and (0) that the full expression for alf{]i) is greater 
than zero as required for stability. 

In some special instances, we can resum the expansion in powers of k; a case we will be 
using below is a'f^{kx = 0,ky = 0, kz)- 

.*2^1/2^ 2^ 1 ~ 



- ^ ~ kl h 



[a + 2nuoQ)~^''^ - (a + 2nuoQ + kl/Am^ 



(15) 



where 

a = 2a + ojq (16) 
measures the distance from the Hc2{T) line. 
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The above results are, however, for the bulk conductivity. In the case with boundaries, 
the calculation should be repeated with the Green's function Q which satisfies the boundary 
conditions imposed on the Ginzburg-Landau equation (^. Let us consider an example with 
a boundary condition at 2; = 0: 

ij{x,y,z = {)-t) = 0. (17) 

One could imagine that such a boundary condition might arise from having a defect like 
a twin-plane at 2; = 0. We admit that this is not a realistic model in any sense for such 
a defect. It is just the simplest example we could devise which enables us to make some 
analytical progress. The response function Q is then given by: 

g{x,y,z]x',y\z']t) = G{x,y,z;x',y',z';t) - G{x,y, -z; x' ,y' , z';t), (18) 

which is just the difference between the bulk Green's function and its image. 

Suppose there is a constant current in the z direction being passed through this boundary, 
and that one would like to calculate the corresponding electric field. Since in this situation 
the current and electric field both point in the z direction, only the zz component of the 
conductivity S'-'^^ in the presence of the boundary is required. Furthermore, the electric 
field varies only in the z direction, hence, the x and y variables can be summed over (or in 
momentum space = 0, ky = 0). Thus we need to evaluate 'E^J']{z, z'). As from now on we 
shall only be considering z components, we shall drop the z indices; for example, E, 
cr^zzCk) cr'^'^^(k), Szzzz S, etc. Furthermore, we will denote the uniform conductivity 
cr*^*^(0) simply as a'^^\ 

The calculation of the conductivity (which is a four-point function) involves a product of 
Green's functions. Consequently, when repeating the steps outlined earlier with Q instead 
of G, one can identify four pieces: 

J:^'\z,z') = a^'\z-z') - a^'\z + z') + c;^'\z,z') - <;^'\-z,z'). (19) 

The first term comes from the product G{z — z')G{z — z') and is, of course, the bulk conduc- 
tivity the Fourier transform of which was calculated above, eq. (|15]). The second 
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term arises from G{z + z')G{z + z') and is the image of the bulk conductivity —a^^\z + z'); 
the negative sign resulting from the derivative with respect to z associated with the current 
operator. The third term originates from the cross term G{z — z')G{z + z') and takes on 
the form: 

q(z, z') = — ■ / dr T dv V du 



327r^mc£^ Jg Jo J-i sinh(a;o''") (1 — M^f^)^/^ 



X 



2mrUV ml , ,2 "^c i\2 



Z-Z) + — t-t(Z + Z 



r(l — u^f ^) r2(l+Mf)2 t'^{1 — uvY 



xexp - A z-z'f - (. + z'f . (20) 

[ 1t[1 + uv) 1t[\ — uv) J 

And finally, the fourth term is the image of the third, again with a negative sign. 

The first piece, the bulk conductivity a'^^\z — z'), decays away from the source at z'; the 

second piece, its image, decays away from the image of the source at —z'. The cross term 

<;{z,z'), on the other hand, decays away from both source and image. As a consequence, 

it could only produce a significant contribution near the boundary where it can be close to 

both source and image. However, at the boundary, its contribution is killed off by its image 

—q{—z, z'). Thus, one expects the cross terms to play a lesser role. For this reason, though 

mainly for reasons of calculational convenience, we will drop the cross terms. (Of course, 

the cross terms while cumbersome could be handled in numerical work if quantitative results 

were required.) 

Thus far the calculation has only concerned the superconducting contribution to the 
conductivity. Before proceeding to determine the electric field, the normal component must 
be included. Here we will simply add a normal contribution o"^"-*, neglecting any interference 
between the two. On the scales of interest, the normal conductivity is local. Hence our 
final expression for the zz component of combined conductivity is: 

E{z, z') = — r dK e-'^'"' ie"^'" - e-'^'") k^"^ + \ . (21) 

Since in the present scenario, we know the current and want the electric field, we need the 
inverse of S(z, z'\ the resistivity p(z, z'\ For an operator of the form (pi]), the inversion is 
readily performed: 
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(22) 



(this ease of inversion is the reason why the "cross terms" in eq. (|19]) were dropped.) 
Notice that without the inclusion of 0"'-"^ the normal component, this integral would have 
been divergent. The electric field is then obtained by integrating over the source j{z'): 



Eiz) 



dz' p{z,z') j{z'). 



(23) 



For a constant current j{z') = jo, one finds: 

2 Jo 



sm{kzz) 



(24) 



Inserting the form of a^^\kz) provided in eq. (|15D provides the electric field distribution 
shown in Figure 1. At the boundary the electric field takes on the value jo/cr^^^ Since 
we imposed the boundary condition tplz = 0) = 0, the current there must be purely normal 
in composition. The electric field then decays toward its bulk value Jo/[o"*-"'' + a^'^^]. Within 
the bulk, the superconducting channel becomes available and thus a smaller electric field is 
required to produce the same current. In expression (P^, the length scale associated with 
the decay of the field to its bulk value depends not only on static quantities such as a (which 
is related to the c-axis correlation length ^ ~ (rrica)^^^'^ /2) but also on dynamic quantites 
such as the normal conductivity a^'^\ 

What would have resulted from the hydrodynamic approach? With a negative coefficient 
S (as calculated), eq. would give: 



E z) 



Jo 



COS 



151 



(25) 



yielding unphysical oscillatory behavior (see Figure 1). Even assuming S* > as required 
for a stable hydrodynamic theory eq. (p4D would yield: 



E(z) 



a 



Jo 

(n) _|_ ^{s) 



1 — exp 



S 



(26) 



which increases from an initial value of zero at the boundary and saturates at the bulk value, 
in direct contrast to the result above where the field was largest at the boundary. Another 
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fix is required; the appropriate boundary condition (-E'(O) = jo/cr^"'^) must be applied. The 
correct phenomenology can be achieved in the hydrodynamic method by using the correct 
boundary condition and assuming S to be positive. However, with the latter step, one 
forfeits making contact with some underlying theory, which must be seen as a failure of this 
approach. 

The attraction of the hydrodynamic approach is its simplicity. The calculation con- 
sidered above is cumbersome even for a single boundary and with simplifying assumptions. 
The complications will only multiply as one considers geometries more appropriate for mod- 
eling multiterminal transport measurements. We suggest the following ad hoc procedure 
which has the virtue of simplicity, gives reasonable agreement with the full theory, and which 
could be generalized to other geometries. We approximate the conductivity by: 



in) ^(^) + [^W]2 + ^(n) ^ 
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which has the correct small-A; behavior (o"*^"^ + a^^^ + S k^, as derived from Landau- Ginzburg 
theory), has the correct large-fc behavior ( a^""^), and is positive definite {i.e. stable). With 
this form for cr{kz), eq. (^4]) yields: 

1/2 



hAz) = — ^ < 1 H TT^Xp 









Or(") 


S 





(28) 



which has the correct limits: E{0) = jo/cr^'^^ at the boundary and E(oo) = jo/icr^"'^ + a^^^) 
deep inside the bulk. Furthermore, it provides a simple expression for the dependence of 
the length scale associated with the decay of E to its bulk value. Its similarity to the full 
expression can be judged from Figure 1. 

As multiterminal transport measurements represent one of the key experimental probes 
of superconducting materials, calculations within such geometries are important for relating 
theory and experiment. Steps toward this end were initiated within the phenomenological 
hydrodynamic approach [Q; however, the subsequent calculation of negative viscosity coef- 
ficients within Ginzburg-Landau theory invalidates the hydrodynamic approach. Here, 
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we have presented a computation which while more cumbersome than hydrodynamics has 
the advantage that it begins with a more fundamental theory and leads to sensible results. 

M.A.M. thanks A.T. Dorsey for introducing him to the problems posed by negative 
viscosity coefficients. T.B. thanks A.J. Bray for useful discussions. 



10 



REFERENCES 

[1] M.C. Marchetti and D.R. Nelson, Phys. Rev. B42, 9938 (1990). 

[2] M.C. Marchetti and D.R. Nelson, Physica C 174 40 (1991). 

[3] M.C. Marchetti, J. Appl. Phys. 69, 5185 (1991). 

[4] D.A. Huse and S.N. Majumdar, Phys. Rev. Lett. 71, 2473 (1993). 

[5] L. Radzihovsky and E. Prey, Phys. Rev. B 48, 10357 (1993). 

[6] L-W. Chen and M.C. Marchetti, Phys. Rev. B 50, 6382 (1994). 

[7] C-Y. Mou, R. Wortis, A.T. Dorsey and D.A. Huse, preprint (1994). 

[8] H. Safar, P.L. Gammel, D.A. Huse, S.N. Majumdar, L.F. Schneemeyer, D.J. Bishop, D. 
Lopez, G. Nieva and F. de la Cruz, Phys. Rev. Lett. 72, 1272 (1994). 

[9] G.J. Ruggeri and D.J. Thouless, J. Phys. F 6, 2063 (1976). 

[10] H.R. Mikeska and H. Schmidt, Z. Phys. 230, 239 (1970). 

[11] R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals, 

(McGraw-Hill, New York 1965). 



11 



FIGURES 



;. 1 The solid line shows the electric field (normalized by jo/c"^"^) as a function of distance 
from the boundary (measured in units of (mca;)^^/^/2, the mean-field c-axis correlation 
length) as calculated from eq. (|2^ ) with o"*-*-'(fcz) given by eq. (|I5]) with a^'^^/a^'^^ = 1 



and LUo/a = 2. The dashed line shows E{z) resulting from eq. (|28|). The third line 
shows is the oscillatory "hydrodynamic" result using the same parameters with S < 
(as calculated). 
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